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Abstract 

Variational perturbation theory is used to determine the decay 
rates of metastable states across a cúbic barrier of arbitrary heigl·it. 
For high barriers, a variational resummation procedure is applied to 
the complex energy eigenvalues obtained from a WKB expansion; for 
low barriers, the variational resummation procedure converts the non- 
Borel-summable Rayleigh-Schròdinger expansion into an exponenti- 
ally fast convergent one. The results in the two regimes match and 
yield very accurate imaginary parts of the energy eigenvalues. This is 
demonstrated by comparison with the complex eigenvalues from Solu- 
tions of the Schròdinger equation via the complex-coordinate rotation 
method. 
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Introduction 



A particle in a one-dimensional anharmonic oscillator with a cúbic potential 
is a popular model for tunneling processes with negligible backscattering. It 
has been widely used in the calculation of non-zero temperature quantum 
transition rates in the presence of friction [1, 2]. All such studies have been 
restricted to the regime of high barriers, where saddle-point and instanton 
methods are applicable. A typical expansion around an infinitely high barrier 
for a classical transition rate in a dissipative medium is given in [3]. 

In the last year, the situation has greatly been improved with the deve- 
lopment of a variational perturbation expansion (VPE) for tunneling rates 
[4, 5]. On the one hand, the available tunneling calculations at high barriers 
were extended to lower barriers. On the other hand, the regime of very low 
barriers became accessible, where metastable states become short-lived and 
the decay proceeds by sliding rather than tunneling. The results converge 
rapidly for all barrier heights, as has been demonstrated for the anharmonic 
oscillator with an x 4 -potential [5, 6] 

The purpose of the present paper is to apply the variational method to the 
cúbic potential and to give imaginary parts for the energy eigenvalues at all 
barrier heights, i.e., for all vàlues of the coupling constants A. In this way we 
go beyond earlier results by Drummond [7] who applied Borel-like methods 
to high-order expansions and was restricted to the regime of intermediate 
barrier heights. 

Our inputs for the variational treatment WKB-expansion for a high 

barrier and the Rayleigh-Schròdinger perturbation series when the barrier is 
low. 

We have reformulated the WKB-method in a way that makes it acces- 
sible to computer àlgebra. This permitted us to determine the expansion 
coefficients around the WKB results up to 7 th order in A 2 . These coeficients 
by themselves do not help to improve the results since the semiclassical ex- 
pansion is divergent. The variational approach resums the expansion and 
leads, even at the level of a first-order WKB-correction, to a remarkable 
improvement of the results. 

In the low-barrier sliding regime, i.e., for not too small vàlues of the 
coupling constant A, a first-order VPE yields the imaginary parts of all ei- 
genvalues E n (X) within ~ 10% of the exact result (see Figure 1). This is all 
the more remarkable as no first order in A appears in the bare perturbation 
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series which contains only even orders in A, for symmetry reasons. To seventh 
order in A, the results deviate by ~ 5x 10~ 5 (see Table 1). In the limit of 
infinitely large A, a simple scaling argument shows that the eigenvalues grow 
like kA 2 / 5 . The VPE reproduces this behavior and yields excel·lent results 
for k (see Table 2). 

Let us remark that the cúbic interaction has also attracted a certain 
amount of interest in quantum field theory, since a </> 3 -theory in six dimensi- 
ons shares certain properties with quantum chromodynamics such as renor- 
malizability and asymptotic freedom [8]. In the field theoretical context, the 
instability of the potential due to tunneling is generally discarded by defining 
the theory to be identical to its perturbative expansion, where all results are 
real numbers. The present method should eventually be generalized to field 
theory and help defining the </> 3 -theory beyond perturbation theory. 

1 General Description of Variational Per- 
turbation Expansion 

We first describe the construction of the VPE in the regime of low barriers. 
Recali standard perturbation theory, which is based on splitting H in two 
parts: 

Ho + -f/int , 

whose first term can be analyzed exactly whereas the second is treated as a 
small correction. For an anharmonic oscillator with a quartic potential, the 
conventional splitting into 

H- = + 
2 dx 2 2 ' 

hïT = g*\ (i) 

yields, for each eigenvalue, an asymptotic series in g which is numerically 
useless unless g is extremely small. In the VPE the splitting is 

rrVPE _ 1 d2 , ^ 2 2 

U ° ~ 9^2 + o X 
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rrVPE U iL 2 : 4 
#int = o X + X ■ 
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The free part of the Hamiltonian H contains a trial frequency to be optimized 
systematically. The perturbation expansion in powers of H^ E for the nth 
excited energy eigenvalues E^ yields an O-dependent expressions: 

E^ YPE = E^\g,ü,u},N), (3) 

where N indicates the order in g up to which the expansion has been per- 
formed. Since the final result must be independent of 0, an optimal energy 
is found by choosing that value of for which the result has the weakest 
O-dependence. 1 This criterion can be implemented in a number of ways: 
Usually one simply looks for a local extremum of E^ n \g,Çl,u, N), i.e., one 
solves 

-A-E^(g,n^,N) = (4) 

for which defines to be a function of g, N and n. The vàlues 

E^(g,ü(9,N,n),u,N) 

are good approximations to the exact eigenvalues even for N = 1, in parti- 
cular, the exact large-g behavior is qualitatively reproduced (E^ = k g 1 / 3 , 
as g — > oo), and the constant k is astonishingly well approximated (k^T] 3 = 
0.6814, compared with the exact value 0.667986). 

For the quartic anharmonic oscillator, the criterion (4) can be applied only 
to odd orders N. Otherwise no zeros are found. In this case, the minimal 
sensitivity on is given by 

d2 EW(g,n,u>,N) = (5) 



dtt 



For large iV, both conditions have many solutions. Then an optimal 
approximation is obtained by searching for plateaux in 0, where the 0- 
dependence is very weak. These plateaux may contain very smooth local 
extrema [6]. If the extrema are plotted in the N- 0-plane, the lower envelope 
of all extrema gives freqüències with the fastest convergence (exponentially 
fast even for infinite coupling strength). 

In this paper, a similar treatment will be given for the cúbic anharmonic 
oscillator, where the asymmetry of the potential requires a few modifications. 



1 This has been called the Principle of Minimal Sensitivity [9]. There exists also a 
slightly different, closely related principle of Fastest Apparent Convergence [10]. 
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In Section 2.2, the O-dependent energies (3) are determined from the con- 
ventional perturbation coefficients by a simple extension of the reexpansion 
trick developed in [5]. 

A convergence proof for the VPE method applied to Rayleigh-Schròdinger 
(RS) perturbation expansions based on a dispersion relation is given in [11] 
(see also [12], where the convergence of the variationally resummed expansion 
for the partition function is proved). The observed convergence is, however, 
much better than the bound given in the proof (cf. [5]). It is exponentially 
fast even for infinitely strong couplings. 



2 Cúbic Oscillator 

We now consider the anharmonic oscillator with a cúbic potential and a 
Hamiltonian: 

1 d 2 w 2 , , , 

H = — + — x 2 -Xx 3 . 6 

2 dx 2 2 y J 

We choose energy units such that u = 1. Due to the asymmetry of the 
potential, the Hamiltonian must be split in a more general way than in (2): 

H = H + XV, (7) 

with 

1 ^ 02 2 o2 

ün = : — ~ ~\ X +JÍ Xn X 

2 dx 2 2 

V = — - x x + —— x 2 + x d . (8) 
A 2 A 

A second variational parameter x permits us to choose an optimal position 
of the potential minimum. The VPE-improved eigenvalues are evaluated by 
solving the following set of equations for x and 0: 

j^El n \\,íl,x ,N) = 

-^—E (n \\, 0, x , N) = 0. (9) 

O Xq 

and inserting the solutions into 

£W VPE (À) = EW(A, 0(A, N, n), x (\, N, n), N). 
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Excel·lent approximations to the complex eigenvalues oí the system are obtai- 
ned by choosing the proper complex zeros of Eqs. (9). The functions x (X) 
and 1 — 2 turn out to be functions of order O(X) and (9(A 2 ), respectively 
so that V can be regarded as being of order A°. 

As a motivation to this splitting we recali that the VPE can be extended 
from energy eigenvalues to path integrals [5, 13, 14], where the lowest ap- 
proximation reduces to the Feynman-Kleinert variational approach (FKVA) 
[15, 16] which is a powerful tool for the approximate calculation of partition 
functions, particle distributions, etc. The path integral in that approach 
contains an integration parameter, the path average 

1 fP 
x o = T, / x(t)(It, 
p Jo 

to be integrated over at the end. The partition function of a quantum me- 
chanical particle 



Z = J T>x exp — J 



rP í x{ T f 



+ V(x(T))\d' 



(10) 



is rewritten as an integral over x : 

Z = 



~P Veff(^o) 

V2^ 



dx . 



:iv, 



with the effective classical potential V e s(x ) defined by the restricted path 
integral 



1 



-Çi 14ff(^o) 



J T>x 8 [^xq 



1 rP 



/3 Jo 



x[t ) dr 



exp 



exp 



P íi(r) 2 O 2 . . . . 2 ] / 

P í O 2 1 

l-—(x(T)-x ) 2 + V[x(r)]\dT 



r m 



The second exponential is treated as a perturbation, is optimized and a 
highly improved expansion is obtained. 
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The ground state energy E^ of the particle can be calculated from Z via 
the relation 

E (o) = _ lim Í5_£. (13) 

But in this limit, the integral (11) can be evaluated in the saddle-point- 
approximation, which means E^ is simply equal to the minimumof V e s(x ). 
Now, V e s(xo) differs from the classical potential in that it contains informa- 
tion about quantum mechanical fluctuations, which tend to smear out the 
classical potential. This means barrier heights in the effective potential are 
lower than in its classical counterpart 2 . If the classical barrier is low enough, 
it will be even completely removed by the fluctuations, and there is our point: 
As the coupling constant A of the cúbic oscillator is increased, the barrier 
height decreases and at a certain point, the minimum of the well will dis- 
appear from the effective potential. Then, the corresponding saddle-point 
in (11) will move into the complex plane and we shall obtain the desired 
complex eigenvalue. 

Obviously, the effective classical potential of the cúbic oscillator can be 
found from a variational approach only if the initial harmònic oscillator Ha- 
miltonian H has, in addition to the trial frequency 0, a variable bottom 
position as in (8). 

This argument can be generalized to include excited states [13]. 



3 The Sliding Regime 
3.1 The First Order in A 

We perform the Rayleigh-Schròdinger expansion of (7) to seventh order in 
A. For greater transparency, we give first the lowest-order results. The 
associated VPE equations are then particularly simple and the optimal vàlues 
of and x are good initial vàlues when searching for extrema at higher 
orders. 

The first order approximation to E^ is: 

E™(\, 0, * , 1) = (n + 1/2) (l±^ - ^ + | - Axg. (14) 

2 In the case of the hydrogen atom, V e jj loses the 1/r-singularity. This is an appealing 
way to explain from a new perspective the stability of matter ([17], ch. 12). 
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The optimization conditions (9) amount to 



5 - + 36(n + 1/2)A 2 



0, 



1 - y/l - 36A 2 (n + l/2)/Q 



= Xq. 



(15) 
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if we use as a supplementary condition that x should be zero when A = 0. 
The proper solution of Eq. (15) is found by following the evolution of the zero 
which is closest to = 1 for small A. For n = all results are real as long as 
A < 0.18 beyond which value it acquires an imaginary part. This is due to 
the fact that the imaginary part of the energy — which can be interpreted as 
a tunneling rate — vanishes exponentially for À — > 0. Obviously, Eqs. (15), 
being polynomial in nature, cannot reproduce exponentially small terms. 
Instead they drop to zero. We can give a physical interpretation of this 
behavior. The VPE is not capable of describing genuine tunneling, whose 
rate has the behavior e~ c /' A ' and possesses no power series expansion. It can 
only account for the process of sliding down into the potential abyss, which 
sets in when the real part of the eigenvalue under consideration lies above 
the potential barrier. 

Figure 1 shows the ground-state results for perturbation orders N = 
1,3,5,7 at the sliding-tunneling-divide. Just as in the case of the quartic 
perturbation, Eqs. (15) do not yield useful results for even N. The higher 
order results do not necessarily drop to zero for small A, but cease to converge 
to the correct result beyond the sliding regime. In accordance with this, Fig. 
2 shows that our approximation yields good approximations closer to A = 
as the quantum number rises. 

3.2 Higher Orders 

On dimensional grounds, the conventional Rayleigh-Schròdinger series for 
the nth excited state is an expansion in powers of A/cu 5 / 2 : 




(16) 
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All odd coeficients vanish. The coefficients are most easily calculated via 
the following recursion relation [7, 18]: 



U - m) B% = \ 



m — 1 



J-3 



EdW f>( n ) 
n k,2 n m-k, 3 

(17) 



(n) 



n j,2 



k=l 



for states with even quantum number n, and 



(j - m) B 



(n) 

m,j 



1 



m — 1 



k=i 



(n) 



-3^,3 



m — k,j 



:i8) 



for odd n. 

It is easy to show that the reorganized series corresponding to the splitting 
(8) is generated from such a series by the following operation generalizing 
the trick in Ref. [5]: 



g A a 9/9 fi 2 e -\bd/8j 



ofi-l^V /4 E^ 



A 2 



o 4 



J 2v 



i/=0 



^0 5 (1 - 12 j A/O 4 ) 5 / 4 



(19) 



wi 



ith 



o 2 - yo 4 - 12 j a 

6A ' 



j := ü 2 x . (20) 
Upon expansion to the desired order in A, a and b have to be replaced by: 

1 -O 2 



A 



J_ 
A' 



This shows that the only input needed for a VPE calculation are the coeffi- 
cients of the divergent series (16). The VPE performs a resummation of this 
series which turns out to be convergent. 
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4 The Tunneling Regime 



For large barriers, i.e., for small A, the imaginary part of the energy eigenva- 
lues is found as usual from the WKB-wave-functions of the cúbic oscillator, 
matched asymptotically to harmònic oscillator wave functions around the 
origin. For the cúbic potential 

üj 2 

V(x) = — x 2 — X x 3 , 
this well-known method [19] leads to the following result: 



lmE^(X) =: £wkb(A, n) 



8 n co íio 5 V +1/2 



7r n\ \X 2 



exp 



2u 5 
Ï5A2 



(21) 



which is correct to zeroth order in A. For the inverted quartic oscillator, the 
WKB-calculation has been pushed to first order in the coupling constant by 
Bender and Wu [18]. Inspired by the method they present in the Appendix of 
[18], we have carried the WKB-calculation for the cúbic oscillator to seventh 
order in A 2 using computer àlgebra (see Appendix A where we described our 
method which is more efficient than the original version by Bender and Wu). 
To first order, we obtain (see Appendix A.2 for the seventh order) 



ImEH(A) = eW KB(A,n)|l-(^-^n-HÍn 2 ) ^ + (9[A 4 ]| (22) 

In Appendix B we have derived the first-order correction factor within the 
path integral formalism [20]. 

The result is now improved by an extension of the variational method 
developed in Refs. [4, 5]. The logarithm of the imaginary part (22) is written 
down for the Hamiltonian H and evaluated at the optimal vàlues of the 
variational parameters and j obtained by the procedure in the last section. 
The VPE equations (9) yield improved vàlues for the imaginary parts of 
the eigenvalues in the tunneling regime [21]. Appendix B shows that this 
procedure is legitimate for the ground-state. For the excited states it turns 
out to match reasonably well onto the results from the preceding section and 
is validated by success. It also compares very favorably with results for the 
cúbic oscillator obtained by more conventional means to which we now turn 
our attention: 
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A 


£>(0) CCR 


£(0)VPE 


0.18 


0.43386176 + 0.02524252 i 


43373 + 02530 i 

V ' • 1 / / 1 / | V ' • V ' L·J 7 / V 7 ' 


1 


61288846 + 40859267 i 


61285 + 40861 i 

V 7 • V / 1 W 7 | V 7 • 1 V 7 \, V / 1 ' 


100 


3 89396500 + 2 82900663 i 

t • V ' • ' t • ' V / 7 V 7 V 7 | L·J • W • / V7 V7 V / V / / ' 


3 89383 + 2 82916 i 


A 


fi(4) CCR 


£i(4) VPE 


0.07 


3.87967181 + 0. 05672784 1 


3.87705 + 0.05811 1 


0.1 


3.59046484 + 0.87272058 i 


3.59583 + 0.87116 z 


1 


8. 13685894 + 5. 82497954 i 


8.13549 + 5.8239 i 


10 


20.4997005 + 14.8887626 i 


20.4963 + 14.8863 i 


20 


27.051325 + 19.651713 i 


27.0469 + 19.6485 z 


50 


39.03 + 28.35 i 


39.0215 + 28.3501 i 



Table 1: Comparison of VPE with the results of a numerical solution of the 
Schròdinger equation via the complex coordinate rotation method CCR for 
the quantum numbers n = and n = 4. The VPE input was the seventh- 
order Rayleigh- Schròdinger expansion. 

5 Comparison with Schròdinger Theory 

A surprising feature of the VPE for complex eigenvalues is that no need arises 
to discuss boundary conditions. This is in contrast to the numerical solution 
of the Schròdinger equation, where a non-Hermitian radiation condition has 
to be imposed. We solved this equation using the complex coordinate rotation 
(CCR) method [22], applying one of the conventional techniques. We have 
chosen to use the elegant formalism of continued matrix fractions [23, 1]. 
Figure 2 shows the result of a VPE calculation based on a seventh-order 
perturbation expansion in A of (19) for the sliding regime, matched to a 
third-order expansion for the tunneling regime. The systematically improved 
WKB results from Appendix A.2 have not been plotted since they almost 
coincide with the variationally improved ones up to the sliding regime, where 
they become inapplicable. Table (1) shows numbers for the seventh-order 
sliding results and their excellent agreement with the numerical Schròdinger 
calculations. 
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n 


K 





0.617159 + 0.448390 i 


1 


2.245672 + 1.631577 i 


2 


4.036069 + 2.932440 i 


3 


6.038479 + 4.387212 i 


4 


8.160972 + 5.929293 i 



Table 2: For large coupling A, the E {n \\) go like E {n \\) -> kA 2 / 5 . The 
table shows the coefficients n for the lowest five states, obtained from an 
ll th -order VPE calculation. 

6 Conclusion 

We have successfully applied the VPE to an asymmetric system with com- 
plex eigenvalues. As a resummation method, VPE is more powerful than 
traditional techniques such as Borel-resummation, Padé-approximants, etc, 
since 

• it gives good results already at the lowest perturbation order , 

• it is vàlid for arbitrarily large vàlues of the coupling constant. 

The asymmetry of the cúbic potential produces a shift of the minimum of 
the effective classical potential as the coupling constant is increased. The 
same situation will be encountered in field theories with symmetry breaking. 

A convergence-proof and an estimate of the remainder term are still mis- 
sing. Both will probably have to rely on the dispersion relation [5, 10]: 

Re£<")(A)=u, (, ! + I) + Ap£ Im£M( , ) _L_. (23) 

A method to reproduce the results from appendix A on a path integral 
level, thus systematically improving on the harmonically matched semiclas- 
sical path integral calculation presented in [24, 25], remains still to be found. 
This is needed to extend the WKB results to quantum field theory. 

It may prové an instructive intermediate step to try and find a recursion 
relation for the coefficients of the WKB-series which is similarly simple as 
the one given by Zinn-Justin in [26] for the systematic corrections to the 
semiclassical energy-level splitting in the double well potential. 
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An immediate application of the present results will be a calculation of 
the quantum mechanical transition rates at non-zero temperature. One can 
exploit the fact that the VPE readily yields complex energy eigenvalues íor 
arbitrarily large quantum numbers. These can be used to find statistical 
averages oí decay rates. We believe that this may ultimately pave the way 
to a similarly efficient variational method íor the computation oí dissipative 
quantum transition rates. 

A Appendix: WKB-Calculation 
A.l The Procedure 

The small-A-limit íor the eigenvalues of the Schròdinger equation 

+ ^-x 2 -Xx 3 -E 4>(x) = (24) 
2 ax z 2 

is related to the % — > 0-limit of the semiclassical approximation by the res- 
caling 

x — > x/X. 

The Schròdinger equation for the cúbic oscillator becomes: 



1 d 2 {u 2 /2)x 2 - x 3 - X 2 E 



2dx 2 X 4 



Íj(x) = 0. (25) 



The asymptotically correct behavior of the eigenvalues as A — > is obtained 
if one notes that the two turning points around the minimum coalesce in this 
limit, which means that, contrary to textbook-WKB, we must not approxi- 
mate the potential linearly in the neighborhood of these turning points but 
have to approximate it quadratically. From a technical point of view, this fa- 
cilitates our task since it means that, in region A of Fig. 3 we can expand our 
solution around harmònic oscillator wave functions instead of Airy-functions. 
In other words: In region A } the wave functions is sufficiently well approxi- 
mated by the Rayleigh- Schròdinger wave function calculated to the desired 
order in A. 

For clarity's sake, we shall present our procedure for the first order in A 2 . 
The result for the seventh order can be found at the end of this appendix: 

M*) = M*) + A^i(z) + 6^ 2 (x) + 0[A 3 ]. 
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The fact that E scales with A 2 is advantageous in another respect as well: 
In region B we have to calculate the WKB-solution of Eq. (24) to first 
WKB-order. Its exponentially decreasing part is: 



4>wkb{x) = p(x) 



-1/2 



exp 



p(0 



- 1/2 —p(0 



-1/2 



de 



(26) 



wi 



ith 



p(x) = \/ uj 2 x 2 — 2 A x 3 — 2 E. 



and E has to be given correctly to first order in A 2 . From perturbation theory 
we obtain: 

E^(X)=üj (n + l) --^-(H + 30n + 30n 2 ) + (9[A 4 ]. 
V 2/ 8 Lú 

In this form, V'wkb contains unappealing elliptic integrals. The expressions 
can be readily processed by computer àlgebra if the WKB-integrand is ex- 
panded to second order in E before performing the integration. Our rescaling 
shows that this renders the WKB-function with sufficient accuracy. Follo- 
wing Bender and Wu [18] we shall neglect the exponentially increasing part 
of the WKB-wave-function at the eigenvalue. Furthermore, the integration 
constant for the integral in Eq. (26) is chosen in such a way that, beyond 
the far (linear) turning point, the integral becomes purely imaginary up to 
terms that vanish for x — > oo. This solution is not a good approximation to 
the true wave function near the far (linear) turning point, but, at a sufficient 
distance, it is vàlid on both sides since it may be continued analytically from 
left to right along a semi-circle in the complex plane, thus keeping at bay the 
singularity at the turning point (see Fig. 3). 

Away from the potential well and its coalescing turning points — i.e., for 
x ~ k/X see Fig. 3 — iprs(%) an d Vhvkb^) match to second order in A: 



V> rs (£/A)=CVw KB (£/A) + (9[A a 



(27) 



where: 



c 



1 + 



\| (A 2 )^ 1 / 2 ^™! 
449 1553 



96 



144 



341 
"48" 



n 2 + 



41 
36 



Í28) 
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The constant k appears only in the next order in A 2 . Knowing the constant 
C we are done. The imaginary part of the eigenvalue can be calculated from 
the quantum mechanical current: 

\ { (] (] \ I rco 

ímE = lim — i>(x)—i>*(x) - r(x)—^(x) / / \i>(0\ 2 d£ , (29) 
^^oo 4i \ dx dx ) I J-oo 

In the numerator of Eq. (29) we may use the WKB wave-function in the limit 
x — > oo where it becomes exact and in the numerator ip may be replaced by 
ip TS , the error being exponentially small. Hence: 

lmE = C 2 / < VgVVs > +0[X% (30) 

and the result in Eq. (22) follows now to second order in A 2 , because: 

A 2 

< V'rsIV'rs >= 1 + ^^(87 + 256 n + 246 n 2 + 164 n 3 ) + 0[A 4 ]. 
A.2 The Seventh-Order Result 

For simplicity we set u = 1. The seventh-order systematic correction to 
e\VKB(A,n) is then given by 

Im£(")(A) = e W K B (A,n)(l + A ;i A 2 + A ;2 A 4 + A ; 3A 6 + 

+k 4 A 8 + h A 10 + h A 12 + k 7 A 14 + C[A 16 ]), 



where 



gra , y \ 1+1/2 

£\VKB(A,n) = — exp 

'irnl WJ 



15A 2 J 



ari' 



d 



169 201 141 2 

L· = n n 

16 8 8 

44507 23639 81 2 12907 3 19881 4 

ko = n n H n -\ n 

512 128 4 64 128 

_ 86071851 24682861 20228607 2 287481 3 
3 ~ 40960 4096 4096 + 256 + 
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5847253 , 356307 - 934407 fi 

H ri H ri ri 

2048 1024 1024 

189244716209 83086372869 52650975367 2 
2621440 327680 U 163840 

4279463287 3 4119876039 4 842550663 5 

n H n -\ n 

32768 65536 16384 

69299813 6 119027547 7 131751387 8 

n n -\ n 

8192 8192 32768 

128830328039451 261646543133639 423592869141879 
41943040 20971520 U 20971520 

19715332860521 3 15339960598811 4 1401355481021 
1310720 U 5242880 n + 524288 ' 

2560290063173 6 25460831311 7 122374368543 8 

H n n n + 

2621440 65536 524288 



n 



31204521765 „ 18576945567 in 

H n n 

262144 1310720 

1027625748709963623 232890795508816773 
6710886400 335544320 

Unmi24048881887 2 664185578830753997 3 
83886080 n 503316480 n 

214792314439921013 4 1384181879205177 5 
335544320 n 41943040 n + 

6816589779335 fi 341591367423137 7 

H n ò H n 7 - 

65536 20971520 

1077894391687337 8 47307658671163 9 
83886080 n 12582912 n + 

16009695476127 in 6433288475823 n 

+ n 10 n 11 + 

5242880 10485760 
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873116441649 12 
20971520 

933142404651555165943 2302284992253612495297 
107374182400 53687091200 

34667142800220818412209 2 55618044990614950369 3 
375809638400 n 503316480 n 

81735871758067504249 4 13377851713824147883 5 
1073741824 n 536870912 n + 

953124390791341475 fi 64443774954315701 7 

H n ò H n 7 + 

536870912 16777216 

1044173780396492593 8 791622589744759973 9 
+ 9395240960 n 2013265920 n 

7250316957748797 in 1631225395815011 n 

n 10 + n 11 - 

134217728 20971520 

6912202077410247 12 397936750990707 
335544320 167772160 

123109418272509 14 
n 



1174405120 

B Appendix: Path Integral Calculation 
of Corrections to Semiclassical Limit 

The n = 0-result in Eq. (22) is reproduced to first order in A 2 by expanding 
the exponential in (10), with the cúbic oscillator potential inserted, around 
an instanton solution to second order in A. Two correlation functions with 
respect to an instanton Green function (cf. [20]) have to be calculated. Using 
computer àlgebra one obtains up to terms which vanish exponentially as 
fi — > oo: 

i rp/2 pp/2 ii finí 

Tx := - / dt dt' < x(tfx(t'f >=-(3 



2 J-/3/2 J-/3/2 8 336 
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1 1-13/2 1-13/2 KKÚ 

T 2 :=- — dt dt' <x , ^t)x(t)x(t r ) 3 >=- — . (31) 

A c \ J-/3/2 -'-(3/2 56 

Part of the first value belongs to the real part of the partition function and 
has to be subtracted. This contribution is found by calculating the first 
correlation with the Green function corresponding to the trivial classical 
solution x c i ttr (t) = which is: 



G tr (t,t') = — e 

Z Lú 



,\t-f\ 



(32) 



The result is 
In the end we obtain: 



1 i-f3/2 ,-f3/2 H 

- dt dt' < x{tfx{t'f > tr = — f3 

2 J-p/2 J-P/2 8 



29 
24' 



7\ + r 2 - t 3 
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calculating the classical action and the fluctuation determinant (see [5], ch. 
17) one obtains to first order in A 2 : 



Lú \Lú° 



vr v A 2 GXP 



Im£(°'(A) 
in agreement with (22) for n = 0. 



2uj 5 
Ï5A 2 



(l-^A 2 )+(9[A 4 ] (33) 
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£(°)(À)/e W KB(A,0) 




Figure 1: Convergence of the VPE for the imaginary part of the ground state 
energy around the tunneling-sliding divide. The dashed line is a numerical 
result obtained by complex coordinate rotation (see Section 4) 
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Im£M(A)/ewKB(A,n) 




0.2 0.4 0.6 0.8 1.0 1.2 



Figure 2: Imaginary parts of energy eigenvalues for quantum numbers 
n = 0,1,2,3,4 divided by zeroth-order WKB results. Sòlid lines on l.h.s.: 
variationally improved WKB results from a 3 rd -order expansion in A (5 th or- 
der for n = 0). On r.h.s.: improved RS results from a 7 th -order expansion. 
Dashed lines: exact results from complex coordinate rotation. Dotted lines: 
l st -order WKB result without variation (2 nd -order for n = 0). 
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Figure 3: Tunneling Rate for High Barriers by WKB 
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